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ABSTRACT 



A frequency-hopping binary frequency shift keying (BSFK) ratio-statistic receiver 
with multihops per data bit is an effective electronic counter-countermeasures (ECCM) 
system against partial-band jamming interference. Interference is modeled as Gaussian 
noise. Orthogonal binary signaling and independent fading diversity is considered over 
frequency-nonselective, slow fading Rayleigh, Rician, and Gaussian channels. A forward 
error correcting coding scheme is implemented for a 1/2 rate convolutional code 
algorithm. The probability of bit error is examined for different levels of diversity, 
thermal noise, severity of fading, fraction of bandwidth jammed, and jamming power. 
Uncoded and coded system comparisons are done to determine worst case performance. 
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I . INTRODUCTION 



A frequency-hopping (FH) binary frequency-shift keying 
(BFSK) with multihops per bit is an effective electronic 
counter-countermeasures (ECCM) system. It provides low 
probability of intercept (LPI) and antijamming capability (AJ) 
for the communication channel. An effective strategy called 
partial-band noise jamming is used against this communication 
system [Ref . l:pp.47l] . 

The principal requirement of the communication system is 
to transmit the source information of a binary data sequence 
by means of M-ary FSK over the channel. We consider that the 
channel is modeled as a fading dispersive channel due to 
multipath problems and the jammer as an additive Gaussian 
noise. The addition of a frequency-hopping scheme provides 
some immunity against interference over an allowed system 
bandwidth of W Hz . 

Many different receiver structures for multihops per data 
bit of FH/MFSK waveforms transmitted over the channel have 
been considered in order to counter or mitigate the effects of 
partial-band jamming and fading dispersive channels. Examples 
of such receivers are: 1) square-law non-linear combining with 
adaptative gain control, 2) soft-limiter amplitude control, 

3) self-normalizing soft decision, and 4) ratio-statistic 
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combiner . 



Keller and Pursley analyzed a ratio-statistic combining 
receiver in a frequency hopping spread spectrum system in the 
presence of partial-band interference for a channel with 
Rayleigh fading [Ref . 2 :p. 145] . Their paper shows that the 
spread spectrum scheme provides some immunity to jamming and 
interference, and the ratio-statistic combining technique 
improved immunity to partial-band jamming. 

Riley analyzed the performance of a fast frequency hopped 
noncoherent BFSK receiver with ratio-statistic combining over 
a Rician fading channel with partial-band interference 
[Ref.3]. This thesis involves a comparison of the performance 
of envelope detectors versus square-law detectors, and the 
following conclusions were established: 1) ratio-statistic 
combining used in conjunction with diversity limits degrada- 
tion due to fading and partial-band interference and generally 
provides an improvement in overall performance, 2) Rician 
fading and partial-band interference have a significant impact 
on receiver performance, 3) ratio-statistic combining provides 
protection against partial-band interference as diversity is 
increased, and 4) for diversity of three and four, the 
receiver using envelope detection performs better than that of 
a receiver using square-law detection. 

Prior work by the authors [Ref . 2 ] - [Ref . 3 ] analyzed the 
performance of uncoded systems using soft decision receivers 
in the presence of fading channels with thermal noise and 
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partial-band noise jamming. 

The purpose of this thesis is to show the performance of 
a forward error control (FEC) coded system using a ratio- 
statistic combining BFSK receiver with envelope or square-law 
detection in a fading channel with partial-band interference. 
The present work examines the probability of coded bit error 
versus different environment cases such as diversity, thermal 
noise, jamming duty factor, signal to jamming ratio, direct 
signal to noise ratio, diffuse signal to noise, envelope and 
square-law detection. Performance curves are obtained for 
coded communication systems. 



3 



II. COMMUNICATION SYSTEM DESCRIPTION 



A. MODEL. 

The general structure of the communication system to be 
evaluated is illustrated in Fig. 1.1. The system is charac- 
terized by an error probability which depends on the charac- 
teristics of the transmitter, the channel, and the receiver. 
This model is comprised of the following elements: 

• Sequence of binary messages. 

• Encoder. 

• M-ary FSK modulator. 

• Spread spectrum hopper. 

• Fading channel. 

• Partial-band interference modeled as additive Gaussian 
noise . 

• Thermal noise. 

• Spread spectrum dehopper. 

• M-ary FSK demodulator with square-law or envelope 
detection. 

• Ratio-statistic soft decision. 

• Decoder. 

• Error corrected binary sequence message. 
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Figure 1.1 Communication System Model 



B. PARTIAL-BAND INTERFERENCE. 

In an interference environment, frequency hopping spread 
spectrum (FHSS) systems are used to reduce receiver perfor- 
mance degradation due to different type of jammers. In order 
to explain the partial-band interference term, it is necessary 
to understand the barrage noise jammer. This jammer transmits 
bandlimited Gaussian noise with a one sided power spectral 
density (PSD) of Nj W/Hz and the jammer power spectrum covers 
exactly the frequency range of the spread spectrum signal. 
The effect of the barrage noise jammer on the system is to 
increase the Gaussian noise at the output of the receiver. If 
FHSS is present, the jamming power is more efficiently used by 
transmitting all the available power in a limited bandwidth 
which is smaller than the spread spectrum signal bandwidth 
[Ref . 4 :pp . 558 ] . This strategy is called partial-band inter- 
ference and the fraction of the spread spectrum signal band- 
width which is jammed is denoted by 7 . For total jamming 
power J and a FHSS signal bandwidth of W, the barrage noise 
jammer one-sided PSD is Nj = J/W over the entire band, while 
the partial-band jammer one-sided PSD is N } /y = J/ 7 W over a 
bandwidth 7 W. The total two-sided noise power spectral den- 
sity is 7' 1 Nj/2+N 0 /2. The partial-band interference is 
assumed to be present in each branch of MFSK demodulator, and 
the fraction 7 of the spread spectrum corrupted is assumed to 
be the same for an entire symbol. The fraction 7 represents 
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the probability that partial band interference is present in 
all M branches of the receiver, and 1-7 is the probability 
that partial-band interference is not present in all branches 
of the receiver. 



Partial -Band Jamming 











l< — 




>1 





Total Bandwidth 



f 



Figure 1.2 Barrage and Partial-Band Interference 



C. FADING MULTIPATH CHANNEL. 

In the multipath propagation model for the channel, the 
fading phenomenon is a result of the time variation in the 
signal phases. Randomly at times, signals can add destruc- 
tively. When that occurs, the resultant received signal can 
be significantly attenuated. At other times, the incoming 
signals add constructively so that the received signal is 
large. These amplitude variations in the received signal are 
called signal fading and are due to the time variant multipath 
characteristic of the channel [Ref . 5 : pp . 457 ] . 

Signal fading appears is produced by the signal bouncing 
off mountains and man-made structures as well as atmospheric 
and other interactions such as refraction and diffraction 
processes. The resultant received signal is modeled as a 
complex valued Gaussian process and is assumed to be comprised 
of a direct and many indirect components. When the impulse 
channel response is modeled as a zero mean, complex valued 
Gaussian process, its magnitude envelope at any time is 
Rayleigh distributed. In this case the channel is said to be 
a Rayleigh fading channel. 

In the event that there are signal reflectors in the 
medium, the channel impulse response is no longer be modeled 
as having a zero mean. In this case, its magnitude envelope 
has a Rician distribution and the channel is said to be a 
Rician fading channel. 
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The probability density function (pdf) of a Rician faded 
signal is [ Ref . 8 : pp. 108 ] 




(l) 



where a>0, a 2 is the power in the direct component, 2a 2 is the 
power in the diffuse components, I 0 (a) is the zeroth-order 
modified Bessel function of the first kind, and u(a) is the 
unit step function. For deep fading on the channel, the 
direct communication path is blocked, and as a consequence the 
direct signal component power is zero. In that event, the 
received signal is Rayleigh distributed, and its pdf is 
[Ref . 6 :pp. 108 ] 



In this thesis, the channel model is assumed to be:l) slow 
fading which implies the signal amplitude remains constant at 
least for the duration of a single hop, and 2) non-frequency 
selective which implies that the hop bandwidth is small 
compared to the coherence bandwidth. This implies that each 
hop frequency can be considered as if it were from an indepen- 
dent channel where the phase and amplitude are relatively 
constant. 



f A (a) = exp [ - ] u(a) 

n 2 z>rt 2 



O 2 2ct 



( 2 ) 
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D. BINARY RATIO-STATISTIC COMBING RECEIVER. 

Ratio-statistic combining implies that decision statistic 
for the receiver is a ratio. The general structure of the 
binary ratio-statistic combiner to be considered is depicted 
in Figures 1.3 and 1.4 [ Ref . 3 : pp . 2 6-2 8 ] . As illustrated, the 
incoming signal is dehopped and nocoherently demodulated by a 
bandpss filter of bandwidth B=R h (where R h is the hop rate) 
followed by an envelope or square-law detector for each of the 
two BFSK frequencies. The ratio-statistic for each detector 
output of a particular diversity reception is equal to the 
detector output divided by the value of the maximum detector 
output. Consequently, one detector output of the ratio- 
statistic receiver for a given diversity reception is always 
unity, and the other detector output for that diversity 
reception is less than unity [ Ref . 2 : pp . 14 6 ] . 




Figure 1.3 Diagram of Ratio Statistic Combiner. 
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Figure 1.4 Alternate Implementation of the Ratio-Statistic 
combining receiver. 

The mathematical analysis to develop the probability 
density function for f w (w) was performed by Riley [Ref.3:pp.27 
-53] where it is shown that Figure 1.3 is equivalent to Figure 
1.4. Furthermore, it is shown that, in Figure 1.4, the random 
variable X lk and X 2k are both normalized to obtain the random 
variable 



ijc “ or* 



(3) 



where i= 1 or 2 and q k = max(X lk , X*) • The random variable Z* 
varies between 0 and 1. The random variable Y k varies from -1 
to +1, and the random variable W varies from -L to +L. 
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The probability of bit error P,^ for a diversity L and a 



jamming duty factor 7 is 




( 4 ) 



where P be (j) represents the conditional probability of bit 
error given that j of L hops are jammed. The number of hops 
that are jammed will depend on the duty factor of the jamming 
7. The total probability of bit error is found by summing 
over all possible values of j. P be (j) is determined by 



The derivation of the PDF for Y k and W requires conside- 
ration of two cases: when the partial-band interference is 
present and when it is not. The evaluation of f w (w) requires 
an L fold convolution. The noise power at the receiver when 
partial-band interference is present with probability 7 is 
given by 




( 5 ) 



o 2 k = (y-'-Nj+Nj B 



( 6 ) 
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The noise power at the receiver when no partial-band 
interference is present with probability 1-7 is 

o 2 k =N 0 B (7) 

The pdf for the random variable Y k is obtained as follows. 

The dehopped signal is expressed as 

S( t) =i/2acos (to t+0) t®) 

where co is the frequency of the dehopped signal, 0 is the 
phase, (2) 5 a is the amplitude, and 0<t<T b 

If the signal is assumed to be in the upper branch of the 
receiver, then f xik is Rician ditributed and is conditional on 
a being present. The pdf f x2k of X 2k in the lower branch is 
Rayleigh distributed. 

The random variables X* are independent, and the joint 
density function f X ikX 2 k( x ik x 2 kl a ) i- s obtained as the product of 
f xlk and f X2 k • Next, appropiate auxiliary variables are 

defined, and the Jacobian of the transformation is obtained. 
Finally, the conditional pdf f Yk (Y k |a) of Y k is obtained by 
integrating over the auxiliary variable. 
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The unconditional pdf for Y k is 



■Yk 



< y *> *f 



f Yk ( Y k \a) f A (a) da 



where f A (a) represents the pdf of the signal amplitude, 
is a Rician random variable as given in equation (2) 
diffuse signal-to-noise ratio is defined as 



The direct signal to noise ratio is defined as 




and the direct to diffuse ratio as 



DDR= 




The PDF for the random variable W as f w (w) 
evaluated from 



is numer 



fjm=[fy i( Y k ) ] 0J '® [ fy° k ( y k ) ] 01 [L ' j) 



(9) 

which 

The 

( 10 ) 



( 11 ) 



( 12 ) 

ically 

(13) 
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where the j superscript implies jammed hops and the o super- 
script implies un jammed hops. 

P^j) is obtained by equation (5), and the probability of 
bit error is found by numerically evaluating equation (4). 



1. Unconditional pdf fvk(Y k ) for envelope detection 

The conditional pdf for the random variable X lk is 

given by 



UiJa) = 



-Xlk V ^1 k 



IT 



(X? k +2a 2 ) 



20* 






v/2 aX 



lJc 



( 14 ) 



The pdf for the variable X 2k is 



(X, J =. x * k 



-X2k'^2k 



20 * 



ol 



( 15 ) 



The probability density function for the random 
variable Y k in the interval of -1 < Y k < 0 is given by 
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fyk ( ^ ^ ~ ' 



2 (i + yj 



Pk 

i+i k +(i *Y k ) 2 



X 



1 + 



[1+ (1 + Y*) 2 ] d+5*+(l + r*) 2 ] 

(1+Y k ) 2 ( p^d+d+r*) 2 ) ^ 



( 16 ) 



i + S* + d+i *) 2 



5i 



i+d+d+Yj 2 



/j 



and for 0 < Y k < 1 



fy k (y K )=- 



2 d-yj 



d+d-Yj 2 ] [i + (5v+D d-YJ 2 ] 



PlrU -^) 2 
i*(? k +D (i -y k )'- 



(17) 



1 + - 



1+dv+l) ( 1 - Y k ) 2 






p^(H-d-Y^) 2 ) 
l + d*+l) (1-Y*) 2 



/J 



2. Unconditional pdf f Yk (Y k ) for square law detection. 

The conditional pdf for the random variable X lk is 

given by 



Ui k +2a 2 ) 

f x M k l a )=A e 2 ° l J o ( 

ol 



</2aX lk) 



(18) 



O ' 5 

o 2 



The pdf for the random variable X 2k is 



xlk 



W*a*>=-T e 2 °" 



(19) 
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The probability density functions for the random variable Y k 
in the interval 1 < Y k < 0 is given by 



f yk (Y k ) 



x 



1 + 



1 e 2 *( k *y k 

(2 + Y k ) (2 +Z k +Y k ) 

( i , p k^+Y k ) y 
2+Cjk+M * 2 +i k +Y k ) 



( 20 ) 



and for 0 < Y k < 1 



fy k (YJ 



x 



1 + 



. 1 e 2-Y k *l k U~Y k ) 

(2 -Y k ) [2-Y k +t k (l-Y k ) ] 

1 ( , , P k ^-Y k ) V 

2-^ + «jt<l-^) { k 2-Y k +t k {l-Y k ) ) 



( 21 ) 



E. FORWARD ERROR CORRECTION CODING. 

Forward error correction (FEC) is a method that employs 
the adding of systematic redundancy and noise averaging at the 
transmit end of a link such that errors caused by the channel 
can be corrected at the receiver by means of a decoding 
algorithm. The binary data source generates information bits 
at R s seconds. These information bits are encoded at a coded 
rate R, and the coded bit energy is 



E C =R E b 



( 22 ) 
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The encoder output sequence is modulated and transmitted over 
the communication channel. At the receiver, the demodulator 
output is passed through a decoder to recover the original 
binary data. 

FEC can be implemented with either block or convolutional 
codes [Ref . 7 :pp. 417-441] . Convolutional codes are selected 
for this thesis. 

A. J. Viterbi [Ref. 8] defines a convolutional encoder as 
a linear finite state machine consisting of a K-stage 
shift register and n linear algebraic function generators. 
The binary input data is shifted along the register b bits at 
a time. 

Each coded symbol carries an average of b/n information 
bits, and the code rate R is 

R=— (23) 

n 

Figure 1.5 depicts the general structure of a convolu- 
tional encoder, where for a given information bit rate the 
required E b /N 0 for a specific bit error rate with FEC is 
generally less than the E b /N c required without FEC. Coding 
gain is defined as the difference in E b /N 0 with and without 
FEC [Ref .7:pp.419] . 
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CONSTRAINT LENGTH K-7 




COOED OUTPUT DATA 



Figure 1.5 A Convolutional Encoder. 

The Viterbi decoder is one of the most common decoders for 
convolutional codes because it is the optimum decoding 
algorithm in the sense of maximum likelihood decoding of the 
entire sequence for convolutional codes. The Viterbi decoder 
performance can be approximated through the use of an upper 
bound [Ref . 9 : pp . 407 ] . For the rate of 1/2, constraint length 
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K=7 convolutional code, the decoded bit error probability is 
bounded by 



p dbe < \ (36£ 10 + 211D 12 + 1404D 16 +11633D 18 + ) ( 24 ) 

where 

d = 2 /p^n^r (25) 
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III. METHOD OF RESEARCH 



Computation of the probability of bit error involves a 
numerical evaluation of (12), (33), and (34). Numerical 
evaluations are done for both envelope and square-law 
detection. 

The coded system performance is shown by plotting the 
decoded bit error probability P dbe versus the average signal 
bit energy density-to- jamming noise power spectral density 
ratio E b /Nj. Several constraints are established in order to 
analyze the performance of the communication system. 

The jamming duty factor y takes the values .001, .05, .25, 
and 1.0 to determine the effects of different amounts of 
partial-band interference on the system. The direct-to- 
diffuse ratio DDR takes the values .01, 10, and 10000 to 
determine the effects of Rayleigh, Rician, and Gaussian 
channels on the system. The bit energy-to-noise ratio E b /N 0 
takes the values 13.35 and 16.00 dB to evaluate the effects of 
thermal noise on the system. The bit energy-to- jammer density 
ratio E b /Nj takes values from 0 to 40 dB to review the effects 
of jamming against the system and each constraint is evaluated 
for diversity L values from 1 to 4 . 

The uncoded and coded performances are compared in order 
to determine the advantages of using FEC with the spread 
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spectrum frequency hopping/noncoherent binary FSK ratio- 
statistic combining receiver in a communication channel 
characterized by fading and partial-band interference. 
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IV. DISCUSSION OP RESULTS 



The channel is modeled for a moderately strong direct-to- 
diffuse signal ratio a 2 k /2a 2 =10 which is considered Rician, for 
a 2 k /2cr 2 =.01 which is considered a Rayleigh faded channel, and 
Gaussian with a 2 k /2cr 2 =10000 . The worst jamming case corres- 
ponds to broadband jamming 7 = 1 . 

Performance with FEC in all cases showed an appreciable 
improvement in performance as expected. Comparison of the 
rate 1/2 convolutional code versus uncoded performance shows 
the expected coding gain. The coding gain is based on an 
optimum diversity level, this being the level where the 
probability bit error approches an asymptotic lower limit for 
a specific value of bit energy-to-noise density ratio. 
Uncoded performance showed an asymptotic lower limit of 
P^IO 4 for a bit energy-to-noise of 16 dB and P^IO ' 3 for 13.35 
dB; however, decoded performances show an asymptotic lower 
limit of P be = 10' 9 for a bit energy-to-noise density ratio of 16 
db and P^IO ' 6 for 13.35 dB. The worst case occurs at the point 
where the probability of the bit error reaches its maximum. 

Figures 3.1 to 3.4 are the envelope receiver performance 
for partial-band and broadband jamming at a direct-to-dif fuse 
ratio of 10 for diversities L=l, 2 , 3 and 4. Figure 3.1 shows 
that for bit energy-to- jamming density ratios greater than 10 
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dB, the performance with broadband jamming is better than 25% 
partial-band interference. For bit energy-to- jamming ratios 
greater than 15 dB, again, the performance with broadband 
jamming is better than both 5% and 25% partial-band inter- 
ference. The uncoded system shows performance degradation due 
partial-band interference from 8 to 40 dB, while the coded 
system shows better performance against partial-band inter- 
ference from 15 to 40 dB. For bit energy-to- jamming density 
ratios less than 15 dB, receiver performance improves more for 
partial-band than broadband jamming. 

Figures 3.2, 3.3, and 3.4 show that degradation due to 

partial-band jamming is reduced dramatically for L=2 , 3, and 
4. Figure 3.5 shows the worst case for diversities of L=l, 3, 
and 4. It is interesting to note that if we compare the plots 
for diversities of L=1 and 4, the better performance 
corresponds to diversity L=4 , but between L=3 and 4, diversity 
L=3 performs better than L=4 . There are improvements in 
performance when diversity L=3 is used for bit energy-to- 
jamming density ratios above 18 dB. 

Figure 3.6 is the receiver performance in the presence 
of a Gaussian channel which is degraded significantly for no 
diversity while Figure 3.7 shows how this degradation is 
reduced when diversity is increased. 

Figure 3.8, for a Gaussian channel, we see that the 
performance advantage due to increased diversity is lost for 
bit energy-to- jamming density ratios greater than 26 dB, where 
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the L=1 worst case performance drops below the L=4 worst case 
performance. 

Figures 3.9 and 3.10 show the receiver performance for a 
Rayleigh channel; the worst case performance corresponds to 
broadband jamming and the effects of increasing diversity are 
an important consideration. 

Figure 3.11 shows, for different amounts of channel 
fading, the receiver performance with 5% partial-band inter- 
ference and diversity L=4 . Figure 3.12 illustrates the 

effects of broadband jamming. The worst case performance for 
Rayleigh channels corresponds to broadband jamming. 

For E b /N 0 =13 . 35 dB, Figures 3.13, 3.14, 3.15, and 3.16 show 
receiver performance for different fractions of partial-band 
and broadband jamming in the presence of a Rician faded 
channel for diversities of L=l, 2, 3 and 4. For no diversity, 
the performance with broadband jamming is better than 5% 
partial-band interference for bit energy-to- jamming density 
ratios from 17 to 30 dB and better than 25% partial-band 
interference for bit energy-to- jamming density ratios lower 
than 30 dB. Partial-band interference degrades the receiver 
performance dramatically versus the receiver performance for 
bit energy-to-noise of 16 dB. When diversity is increased, 
degradation is reduced for diversity L=2 and worst case 
performance occurs for broadband jamming; however, diversity 
greater than L=2 does not improve performance because 
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diversity L=2 has an asymptotic lower limit of P^IO' 5 while 
for L=4 the asymptotic lower limit is P^lO^*. 

Figure 3.17 shows worst case performance for diversities 
L=l, 3, and 4. The advantage of using diversity is lost for 
low values of bit energy-to-noise density ratios. 

Figure 3.18 shows receiver performance for a Gaussian 
channel, which is better than that for a Rician faded channel 
(cf . , Figure 3.16) . 

For Gaussian channels, Figure 3.19 shows that performance 
with increased diversity is dramatically poorer for all bit 
energy-to- jamming density ratios. 

Figure 3.20 shows both bit energy-to-noise density ratios 
for Rician faded channels for 13.35 and 16 dB. For larger 
values of energy-to- jamming density ratios, receiver perfor- 
mance tends to be controled by thermal noise. 

Receiver performance for L=3 is better than L=1 for a bit 
energy-to-noise density ratio of 16 dB for all bit energy-to- 
jamming density ratios above 17 dB. For a bit energy-to-noise 
density ratio of 13.35 dB, lower orders of diversity perform 
better than higher for all values of bit energy-to- jamming 
density ratios. 

In Figures 3.21, 3.22, and 3.23, the performance of the 

square-law receiver is illustrated for a bit energy-to noise 
ratio of 16 dB and diversities of L=l, 3 and 4 for a Rician 
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channel. For diversities greater than L=l, the worst case 
performance tends to be when the interference is broadband. 
For the envelope detection case no performance improvement is 
obtained for diversities greater than L=3 . Figure 3.24 and 
Figure 3.25 show worst case performance for diversities of 
L=l, 3, and 4 for bit energy-to-noise density ratios of 16 dB 
and 13.35 dB, respectively. 

Figures 3.26 and 3.27 show comparisons between both coded 
detections examined for bit energy-to-noise density ratios of 
13.35 and 16 dB. For diversities of L=1 and L=3 , the envelope 
detector performs better than the square-law detector for all 
values of bit energy-to- jamming density ratios. For uncoded 
detection, the envelope detector performs better than square- 
law detector for diversities of L=3 or 4. 
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Probability of Bit Error 



ENVELOPE DETECTOR Eb^No= 16 dB L= 1 RATIO= 10 




JAMMING RATIO GAMMA 
*-.001 
o - .05 
- - .25 
-- 1.00 

Figure 3.1. Performance of a FFH/NCBFSK/RS Receiver with 
Envelope detection, where £^1^=16 dB, Direct 
to Diffuse RATIO=10, and L=l. 
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Probability of Bit Error 



ENVELOPE DETECTOR Eb/No=16 dB L=2 RATIO= 10 




*-.001 



o - .05 
.25 

-- 1.00 

Figure 3.2. Performance of a FFH/NCBFSK/RS Receiver with 
Envelope detection, where Ej/N^ie dB, Direct 
to Diffuse RATIO=10, and L=2. 
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Probability of Bit Error 



ENVELOPE DETECTOR Eb/No=16 dB L=3 RATIO= 10 




# -.001 



o - .05 
.25 

- = 1 .00 

Figure 3.3. Performance of a FFH/NCBFSK/RS Receiver with 
Envelope detection, where Ejj/N^l.6 dB, Direct 
to Diffuse RATIO=10, and L=3. 
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Probability of Bit Error 



ENVELOPE DETECTOR Eb/No=16 dB L=4 RATIO= 10 




* = .001 
o - .05 
-.= .25 
- = 1.00 

Figure 3.4. Performance of a FFH/NCBFSK/RS Receiver with 
Envelope detection/ where £^1^=16 dB/ Direct 
to Diffuse RATIO=10/ and L=4. 
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Probability of Bit Error 



ENVELOPE DETECTION Eb/No= 16 dB RATIO= 10 




DIVERSITY 

* - LI o - L3 - - L4 



Figure 3.5. Performance of a FFH/NCBFSK/RS Receiver with 
Envelope detection/ where E^N^ie dB, Direct 
to Diffuse RATIO=lO, and the worst case for 
L=l, 3 and 4. 
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Probability of Bit Error 



ENVELOPE DETECTOR Eb/No=16 dB L=1 RATIO= 10000 




*-.001 
o - .05 
.25 

-- 1.00 

Figure 3.6. Performance of a FFH/NCBFSK/RS Receiver with 
Envelope detection, where dB, Direct 

to Diffuse RAT!0=10000, and L=l. 
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Probability of Bit Error 



ENVELOPE DETECTOR Eb/No=16 dB L=4 RATIO= 10000 




JAMMING RATIO GAMMA 



*-.001 
o - .05 
.25 

-- 1.00 

Figure 3.7. Performance of a FFH/NCBFSK/RS Receiver with 
Envelope detection/ where E^N^IS dB/ Direct 
to Diffuse RAT 10= 10000/ and L=4. 
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Probability of Bit Error 



ENVELOPE DETECTOR Eb/No= 16 dB WORST L= 1,4 RATIO= 10000 




Eb/Nj 



DIVERSITY 

- LI - - L4 



Figure 3.8. Performance of a FFH/NCBFSK/RS Receiver with 
Envelope detection, where E^N^ie dB, Direct 
to Diffuse RAT!0=10000, and Worst L=1 and 4 . 
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Probability of Bit Error 



ENVELOPE DETECTOR Eb/No=16 dB GAMMA = .05 RATIO= .01 




DIVERSITY 

* - LI o = L 2 - L3 - - L4 

Figure 3.9. Performance of a FFH/NCBFSK/RS Receiver with 
Envelope detection, where E^N^ie dB, Direct to 
Diffuse RATIOS. 01, GAMMAS. 05, and Lsl,2,3 and 4. 
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Probability of Bit Error 



ENVELOPE DETECTOR Eb/No = 16 dB GANHA= 1 RAT 10= .01 




DIVERSITY 

* - LI o - L2 - L3 - - L4 

Figure 3.10. Performance of a FFH/NCBFSK/RS Receiver with 
Envelope detection, where dB, Direct to 

Diffuse RATIO=.01, GAMHA=1, and L=l, 2 , 3 and 4. 
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Probability of Bit Error 



ENVELOPE DETECTOR Eb/No= 16 dB GAMMA = .05 L= 4 




*-.01 0 - 10 - - 10000 



Figure 3.11. Performance of a FFH/NCBFSK/RS Receiver with 
Envelope detection, where E^N^ie dB, Different 
Amounts of Fading, GAMMA=.05, and L=4. 



38 



Probability of Bit Error 



ENVELOPE DETECTOR Eb/No= 16 dB GANNA= 1 L= 4 




DIRECT-TO-DIFFUSE RATIO 
* - .01 o- 10 - - 10000 



Figure 3.12. Performance of a FFH/NCBFSK/RS Receiver with 
Envelope detection, where £^1^=16 dB, Different 
Amounts of Fading, GAMHA=1, and L=4 . 
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m i n i i » — ni r m 



Probability of Bit Error 



ENVELOPE DETECTOR Eb/No= 13.35 dB L= 1 RATIO= 10 




* = .001 



o = .05 
-.= .25 
- = 1.00 

Figure 3.13. Performance of a FFH/NCBFSK/RS Receiver with 
Envelope detection, where £^1^=13.35 dB, Direct 
to Diffuse RATIO=10, and L=l. 
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Probability of Bit Error 



ENVELOPE DETECTOR Eb/No= 13.35 dB L= 2 RATIO= 10 




*-.001 



o - .05 
-.-.25 
- - 1 .00 

Figure 3.14. Performance of a FFH/NCBFSK/RS Receiver with 
Envelope detection/ where £^^=13.35 dB/ Direct 
to Diffuse RATIO=10/ and L=2. 
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Probability of Bit Error 



ENVELOPE DETECTOR Eb/No= 13.35 dB 



L= 3 RAT 10= 10 




*-.001 
o - .05 
.25 

- = 1.00 

Figure 3.15. Performance of a FFH/NCBFSK/RS Receiver with 
Envelope detection, where £^1^=13.35 dB, Direct 
to Diffuse RATIO=10, and L=3. 
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Probability of Bit Error 



ENVELOPE DETECTOR Eb/No= 13.35 dB 



L= 4 RAT 10= 10 




*-.001 
o - .05 
.25 

-- 1.00 

Figure 3.16. Performance of a FFH/NCBFSK/RS Receiver with 
Envelope detection, where E^N^IS.35 dB, Direct 
to Diffuse RATIO=10, and L=4. 
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Probability of Bit Error 



ENVELOPE DETECTION Eb/No= 13.35 dB RATIO= 10 




DIVERSITY 



* - LI o = L3 - - L4 



Figure 3.17. Performance of a FFH/NCBFSK/RS Receiver with 
Envelope detection, where £^1^=13.35 dB, Direct 
to Diffuse RATIO=lO, and the worst case of L=l, 
3 and 4 . 
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Probability of Bit Error 



ENVELOPE DETECTOR Eb/No= 13.35 dB L= 4 RATIO= 10000 




* = .001 
o = .05 
-.= .25 
- = 1.00 

Figure 3.18. Performance of a FFH/NCBFSK/RS Receiver with 
Envelope detection, where £^^=13.35 dB, Direct 
to Diffuse RATI 0=100 00 , L=4 and different amounts 
of interference. 



45 



Probability of Bit Error 



ENVELOPE DETECTOR Eb/No= 13.35 dB UORST L= 1,4 RATIO= 10000 




DIVERSITY 



- LI - - L4 



Figure 3.19. Performance of a FFH/NCBFSK/RS Receiver with 
Envelope detection, where £,^*^=13.35 dB, Direct 
to Diffuse RATIO=10000, and the worst case L=1 
and 4 . 
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Probability of Bit Error 



ENVELOPE DETECTION Eb/No= 13.35 AND 16 dB RATIO= 10 




ENERGY-TO- NOISE RATIO 
13.35 dB 

..-LI --L3 

16.00 dB 

* - LI o - L3 

Figure 3.20. Comparison of the Performance of a FFH/NCBFSK/RS 
Receiver with Envelope detection, where Direct 
to Diffuse RATIO=10, £^1^=13.35, and 16 dB and 
worst L=1 and 3. 
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Probability of Bit Error 



SQUARE LAU DETECTOR Eb/No=16 dB L=1 RATIO= 10 




JAMMING RATIO GAMMA 
* = .001 
o - .05 
-.-.25 
- = 1.00 

Figure 3.21. Performance of a FFH/NCBFSK/RS Receiver with 
Square Law detection, where dB, Direct 

to Diffuse RATIO=10, and L=l. 
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Probability of Bit Error 



SQUARE LAU DETECTOR Eb/No= 16 dB L=3 RATIO= 10 




*- .001 
o * .05 
-.= .25 
- = 1.00 

Figure 3.22. Performance of a FFH/NCBFSK/RS Receiver with 
Square Law detection, where E^N^ie dB, Direct 
to Diffuse RATIO=10, and L=3. 
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Probability of Bit Error 



SQUARE LAW DETECTOR Eb/No=16 dB L= 4 RATIO= 10 




JAMMING RATIO GAMMA 
*-.001 
o- .05 
.25 

-- 1.00 

Figure 3.23. Performance of a FFH/NCBFSK/RS Receiver with 
Square Law detection, where E^N^ie dB, Direct 
to Diffuse RATIO=10, and L=4. 
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Probability of Bit Error 



SQUARE LAW DETECTION Eb/No= 16 dB RATIO= 10 




DIVERSITY 



* - LI o - L3 - - L4 



Figure 3.24. Performance of a FFH/NCBFSK/RS Receiver with 
Square Law detection, where E^N^ie dB, Direct 
to Diffuse RATIO=10, and the worst case of L=l, 
3 , and 4 . 
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Probability of Bit Error 



SQUARE LAW DETECTION Eb/No= 13.35 dB RATIO= 10 




DIVERSITY 

* - LI o - L3 - - L4 



Figure 3.25. Performance of a FFH/NCBFSK/RS Receiver with 
Square Law detection/ where £^1^=13.35 dB, 
Direct to Diffuse RATIO=10, and the worst case 
of L=l, 3, and 4. 
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Probability of Bit Error 



BOTH DETECTION Eb/No= 16 dB RATIO= 10 




ENVELOPE DETECTOR SQUARE LAW DETECTOR 



* - LI 
o- L3 



-.-LI 
- =L3 



Figure 3.26. Comparison of the Performance of a FFH/NCBFSK/RS 
Receiver both detection/ where E^N^IS dB/ Direct 
to Diffuse RATIO=10 / and the worst case of L=1 
and L=3 . 
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Probability of Bit Error 



BOTH DETECTION Eb/No= 13.35 dB RATIO= 10 




ENVELOPE DETECTOR 


SQUARE LAW DETECTOR 


* - LI 


-.-LI 


o» L3 


- =L3 



Figure 3.27. Comparison of the Performance of a FFH/NCBFSK/RS 
Receiver both detection, where E^N^J.3.35 dB, 
Direct to Diffuse RATIO=10, and the worst case of 
L=1 and L=3. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS. 

This thesis has shown the performance of envelope and 
square-law detection of L-hops/symbol FH/BFSK for worst case 
partial-band jamming. The nonlinear ratio-statistic combining 
technique with forward error correcting coding is an effective 
strategy when combined with diversity to mitigate jamming 
effects. 

Ratio-statistic combining soft decision receivers with 
forward error correcting coding outperforms the uncoded 
receiver for all cases of diversity, noise jamming inter- 
ference, and channel fading. 

The envelope detection optimum diversity level tends to be 
L=3 with a bit energy-to- jamming density ratio of 16 db; this 
corresponds to broadband jamming. 

System performance is dramatically degraded when the 
bit energy-to-noise density ratio is reduced to 13.35 dB. The 
worst case performance for no diversity is better than higher 
orders; however, the worst case performance is better than 
that of the uncoded system under the same conditions. 

The coded performance of square-law detection systems is 
not improved by diversity when jamming is broadband for all 
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channel types. 

Comparision of the performance of envelope detectors and 
square-law detectors show that envelope detection performs 
better than square-law detection for all bit energy-to- jamming 
density ratios considered . 



B . RECOMMENDATIONS . 

The FFHSS/noncoherent BSFK receiver with ratio-statistic 
combining and forward error correcting coding has been 
evaluated for partial-band jamming and for different types of 
faded channels. This analysis should be extended to the 
general M-ary case. In addition, the ratio-statistic system 
should be simulated using some of the available communication 
simulation software packages. 
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APPENDIX 



MATLAB SOURCE CODES 

% The following source codes were developed by 

% JOHN RILEY, 1990 and this codes were modified 

% to implement the 1/2 VITERBI decoder algorithm. 

% ENVMAIN.M Envelope Detection. 

% This is the main source code that generates 

% the decoded Pbe for the system with Ratio 

% Statistic combiner and convolutional VITERBI 

% decoder of rate 1/2. 

% 

r = 1/2; 

RATIO = ' INPUT THE VALUE OF DIRECT TO DIFFUSE RATIO '; 
n= 150; 

Eb = 10 A (X. XX) ; ' INPUT THE VALUE OF BIT ENERGY' 

Ec = r*Eb 

No = 1.0; 

GAMMA =' INPUT THE VALUE OF FRACTION JAMMING ' ; 

L = ' INPUT THE DESIRE SYMBOL DIVERSITY ' ; 

Eh = Ec/L; 

DIFFUSE = Eh/ (RATIO+1) ; 

DIRECT = DIFFUSE*RATIO; 

RHOn = DIRECT; 

XIn = DIFFUSE; 
for xx =1:21 

clear Pb 

EbNj = (xx-1) *2 ; 

Nj = 10 A (-EbNj/10) *Ec; 

Nt = 1 + Nj /GAMMA; 

RHOi = DIRECT/Nt; 

Xli = DIFFUSE/Nt; 
envpbe ; 

PPbe(xx) = Pb; 

X (xx) = 10 * loglO (Eb/Nj ) ; 

end 

save PPbe; 
save X; 
clear 

semi logy (X,PPbe) ; 

title ( 'Eb/No= dB GAMMA= L= RATIO= '); 

xlabel ( 'Eb/Nj ' ) ; 

ylabel ( 'Probability of Bit Error'); 
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% ENVPBE.M for rate code of 1/2 

% This program will determine the probability of error 

% as in equation 5 of Keller's paper. It is the 

% algorithm for the 1/2 VITERBI decoder for envelope 

% detection. 

% 

Pbee =0; 
for j = 0:L 
convfad; 

simper r ; 

Pbc( j+l)=Pbej* (fact (L) / (fact (j ) *f act (L-j ) ) ) *GAMMA" j * ( 1-GAMMA 

) - (L-j ) ; 

end 

Pbee=sum (Pbc) 

d = 2*sqrt ( (1-Pbee) *Pbee) ; 
valor (xx) =d 

qq = 36* (d) "10; 
ww = 211* (d) "12 ; 
ee = 1404* (d) "14; 

Af = 11633* (d) "16; 

Pb= . 5* (qq+ww+ee+Af ) ; 
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% FADNOINT.M Envelope detection 

% This program will sample the probability density 

% functionm of the variable y, it generates the NON 

% JAMMED HOPS sequence and put it into and array f (y) . 

% 

clear f 
for x = l:n 

z = x - 1 ; 
y = -l+z/n; 

C = 2* ( 1+y) / ( ( l+(l+y) A 2) * ( l+( 1+y) ~2+XIn) ) ; 

D = RHOn/ ( 1+ ( 1+y) ^2+XIn) ; 

E = (1+y) *2/ (l+XIn+(l+y) ~2) ; 

K = RHOn* ( 1+ ( 1+y ) ~ 2 ) / ( 1+XI n+ ( 1+y ) ~ 2 ) ; 
f ( x) = C*exp(-D) *(l+E*(XIn+K) ) ; 

end 

for x = l:n+l 

z = x - 1 ; 

Y = z/n; 

C = 2* (1-y) / ( (l+(l-y) ~2) * (l+(l-y) ~2* (1+XIn) ) ) ; 

D = RHOn* (1-y) A 2/ (l+(l-y) *2* (1+XIn) ) ; 

E = 1/ (l+(XIn+l) *(l-y) ~2) ; 

K = RHOn* ( 1+ ( 1-y) ~2) / ( 1+ (XIn+1) * ( 1-y) A 2) ; 
f (x+n) = C*exp(-D) * (1+E* (XIn+K) ) ; 

end 

v = -1 : 2/length (f) : l-2/length(f ) ; 
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oY> <X> <X> <X> <X> 



for x 



end 
for x 



end 
v = - 



FADWINT.N Envelope detection 

This program will sample the probability density 
function of the variable y, it generates the JAMMED 
HOPS sequence and put it into and array f (y) . 

clear f 
= 1 : n 

z = x - 1; 
y = -1+z/n; 

C = 2*(l+y) / ( (l+(l+y) *2) *(l+(l+y) *2+XIi) ) ; 

D = RHOi/ (l+(l+y) A 2+XIi) ; 

E = (1+y) * 2 / (l+XIi+(l+y) "2) ; 

K = RHOi* ( 1+ ( l+y) *2 ) / ( i+xii+ (l+y) *2 ) ; 
f ( x ) = C*exp(-D) * (l+E* (Xli+K) ) ; 

= 1 : n+1 

z = x - 1; 
y = z/n; 

C = 2*(l-y) / ( (l+(l-y) * 2 ) *(l+(l-y) ~2* ( 1+XIi) ) ) ; 

D = RHOi* (1-y) A 2/ (l+(l-y) * 2 * (l+xii) ) ; 

E = 1/ (l+(XIi+l) *(l-y) ~2) ; 

K = RHOi* ( 1+ ( l-y) ~2) / ( 1+ (Xli+l) * ( 1-y) ~2 ) ; 
f (x+n) = C*exp (-D) * ( l+E* (Xli+K) ) ; 

1 : 2/ length (f) : 1-2/ length (f) ; 
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% CONVFAD.M ENVELOPE DETECTION 

% This rutine generates the convolution of the Folded 

% convolutions for the jammed and non jammed HOPS. 

% fw represents the numerical evaluated PROBABILITY 

% DENSITY FUNCTION for the random variable W. 

% 

f adwint 

h = 2/ (length (f) -1) ; 
if j==0 

fwi = [ 1 ] ; 
elseif j==l 

fwi = f; 
elseif j==2 

fwi = conv ( f , f ) ; 
elseif j==3 

fwi = conv ( conv (f , f ) , f ) ; 
elseif j==4 

fwi = conv ( conv ( conv (f , f ) , f ) , f ) ; 
elseif j==5 

fwi = conv (conv ( conv (conv ( f , f) , f) , f) , f) ; 
elseif j==6 

fwi = conv (conv ( conv (conv (conv ( f , f) , f ), f) , f) , f) ; 

else 

fwi = error; 

end 

f adnoint 



if (L- j ) ==0 

fwn = [ 1 ] ; 
elseif ( L- j ) ==1 
fwn = f ; 
elseif ( L- j ) ==2 

fwn = conv (f , f ) ; 
elseif ( L- j ) ==3 

fwn = conv (conv(f , f ) , f ) ; 
elseif (L-j ) ==4 

fwn = conv (conv (conv ( f , f) , f) , f) ; 
elseif (L- j ) ==5 

fwn = conv ( conv ( conv ( conv (f , f ) , f ) , f ) , f ) ; 
elseif (L-j ) ==6 

fwn = conv ( conv ( conv ( conv ( conv (f , f ) , f ) , f ) , f ) , f ) ; 

else 

fwn = error; 

end 

fw = h~ (L-l) *conv (fwi , fwn) ; 
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o\° o\° <*° o\o oV> <A° o\° 



SIMPERR.M 

This program uses Simpson's rule to find the area 
under the array fw for y less than zero. This will 
be the probability of bit error for a given # HOPS 
jammed, Pbe(j) in a receiver using Ratio-Statistic 
Combining. The result is in Pbej . 

areal=0; 

Pbej = 0; 
for s = 2 : (length (fw) -1) /2-1 
if round (s/ 2) == s/2 

areal = areal + 4*fw(s) ; 
else 

areal = areal + 2*fw(s) ; 

end 

end 

areal = areal + fw(l) +fw( (length(fw) -1) /2) ; 

Pbej = areal*h/3; 
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% SQMAIN.M SQUARE LAW DETECTION 

% This is the main source code that generates 
% the decoded Pbe for the system with Ratio 
% Statistic combiner and convolutional VITERBI 

% decoder of rate 1/2. 

% 

r = 1/2; 

RATIO = ' INPUT THE VALUE OF DIRECT TO DIFFUSE RATIO '; 
n= 150; 

Eb = 1 0 A ( x . xx ) ; ' INPUT THE VALUE OF BIT ENERGY ' 

Ec = r*Eb 

No = 1.0; 

GAMMA =' INPUT THE VALUE OF FRACTION JAMMING ' ; 

L = ' INPUT THE DESIRE SYMBOL DIVERSITY ' ; 

Eh = Ec/L; 

DIFFUSE = Eh/ (RATIO+1) ; 

DIRECT = DIFFUSE* RATIO; 

RHOn = DIRECT; 

XIn = DIFFUSE; 
for xx =1:21 

clear Pb 

EbNj = (xx-1) *2 ; 

Nj = 10~ (-EbNj/10) *Ec; 

Nt = 1 + Nj /GAMMA; 

RHOi = DIRECT/Nt; 

Xli = DIFFUSE/Nt; 
sqpbe ; 

PPbe(xx) = Pb; 

X (xx) = 10 * loglO ( Eb/N j ) ; 

end 

save PPbe; 
save X; 
clear 

semi logy (X, PPbe) ; 

title ( 'Eb/No= dB GAMMA= L= RAT 10= •); 

xlabel ( ' Eb/Nj ' ) ; 

ylabel ( 'Probability of Bit Error'); 
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% 

% 

% 

% 



% 



SQCORRE.M for rate code of 1/2 

This program will determine the probability of error 
as in equation 5 of Keller's paper. It is the 
algorithm for the 1/2 VITERBI DECODER for the square 
law detection. 



% 

Pbee =0; 
for j = 0:L 



sqcovfad; 



simper r ; 



Pbc ( j + 1) =Pbe j * ( fact (L) / (fact (j ) *f act (L-j ) ) ) *GAMMA*j *( 1-GAMMA 

) * (L“j) ; 



d = 2*sqrt ( (1-Pbee) *Pbee) ; 
valor (xx) =d 
qq = 36* (d) *10; 
ww = 211* (d) * 12 ; 
ee = 1404 * (d) *14; 

Af = 11633* (d) *16; 

Pb= . 5* (qq+ww+ee+Af ) ; 



end 

Pbee=sum(Pbc) 



64 



<N> <*> oV> oV> oY> 



SQFADNOI . M Square law detection 

This program will sample the probability density 
function of the variable y , it generates NON JAMMED 
HOPS sequence and put it into and array f(y). 

clear f 
for x = l:n 

z = x - 1 ; 
y = -1+z/n; 

C = (1/ ( (2+y) * (XIn+2+y) ) ) ; 

D = ( 1+y) / (XIn+2+y) ; 

E = -RHOn/ (XIn+2+y) ; 

K = RHOn* (2+y) / (XIn+2+y) ; 
f ( x ) = C*exp(E) *(l+D*(XIn+K) ) ; 

end 

for x = l:n+l 

z = x - 1; 

Y = z/n; 

C = 1/ ( (2-y+XIn* (l-y) ) * (2-y) ) ; 

D = RHOn* (l-y) / (2-y+XIn* (l-y) ) ; 

E = 1/ (2-y+XIn* (l-y) ) ; 

K = RHOn* (2-y) / (2-y+XIn* (l-y) ) ; 
f (x+n) = C*exp(-D>* (1+E* (XIn+K) ) ; 

end 

v = -1: 2/length(f ) : 1-2/length (f) ; 
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oV> <*> 



% 

for x 



end 
for x 



end 
v = - 



SQFADWIN.M Square law detection 

This program will sample the probability density 
function of the variable y , it generates the JAMMED 
HOPS sequence and put it into and array f(y). 

clear f 
= 1 : n 

z = x - 1 ; 
y = -1+z/n; 

C = ( 1 / ( ( 2 +y ) * ( XI i+2 +y ) ) ) ; 

D = ( l+y) / (xii+2+y) ; 

E = -RHOi/ (XIi+2+y) ; 

K = RH0i*(2+y) / (XIi+2+y) ; 
f ( x) = C*exp(E) * (1+D* (Xli+K) ) ; 

= l:n+l 

z = x - 1; 
y = z/n; 

c = 1/ ( (2-y+xli* (l-y) ) * (2-y) ) ; 

D = RHOi* (l-y) / (2-y+XIi* (l-y) ) ; 

E = 1/ (2-y+XIi* (l-y) ) ; 

K = RHOi* (2-y) / (2-y+XIi* (l-y) ) ; 
f (x+n) = C*exp ( -D) * ( 1+E* (Xli+K) ) ; 

1:2/ length ( f ) : 1-2 / length ( f ) ; 
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% SQCOVFAD.M SQUARE LAW DETECTION 

% The rutine generates the convolution of the Folded 

% convolutions for jammed and nonjammed HOPS. 

% fw represents the numerical evaluated PROBABILITY 

% DENSITY FUNCTION for the random variable W. 

% 



sqfadwin 

h = 2/ ( length (f ) -1) ; 
if j ==0 

fwi = [ 1 ] ; 
elseif j==l 

fwi = f; 
elseif j==2 

fwi = conv(f , f ) ; 
elseif j==3 

fwi = conv(conv(f , f ) , f ) ; 
elseif j==4 

fwi = conv (conv (conv (f ,f) , f) , f) ; 
elseif j==5 

fwi = conv ( conv ( conv ( conv (f , f ) , f ) , f ) , f ) ; 
elseif j==6 

fwi = conv (conv (conv (conv (conv(f, f) , f) , f) , f) , f) ; 

else 

fwi = error; 

end 

sqf adnoi 
if (L- j ) ==0 



fwn = [ 1 ] ; 



elseif 


(L-j) 




fwn 


elseif 


(L-j) 




fwn 


elseif 


(L-j) 




fwn 


elseif 


(L-j) 




fwn 


elseif 


(L-j) 




fwn 


elseif 


(L-j) 




fwn 


else 






fwn 


end 




fw = h' 


(L-l) 



) ,f) ; 

,f) ,f) ,f) ,f) ; 
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